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We introduce an arbitrary number of accessible modes when analyzing bipartite entanglement 
degradation due to Unruh effect between two partners Alice and Rob. Under the single mode 
approximation (SMA) a fermion field only had a few accessible levels due to Pauli exclusion principle, 
conversely to bosonic fields which had an infinite number of excitable levels. This was argued to 
justify entanglement survival in the fermionic case in the SMA infinite acceleration limit. Here we 
relax SMA. Hence, an infinite number of modes are excited as the observer Rob accelerates, even 
for a fermion field. We will prove that, despite this analogy with the bosonic case, entanglement 
loss is limited. We will show that this comes from fermionic statistics through the characteristic 
structure it imposes on the infinite dimensional density matrix for Rob. Surprisingly, the surviving 
entanglement is independent of the specific maximally entangled state chosen, the kind of fermionic 
field analyzed, and the number of accessible modes considered. We shall discuss whether this 
surviving entanglement goes beyond the purely statistical correlations, giving insight concerning 
the black hole information paradox. 

PACS numbers: 03.67.Mn, 03.65.-w, 03.65.Yz, 04.62.+V 



I. INTRODUCTION 

Studying quantum information in non-inertial set- 
tings requires using tools coming from general relativity 
EE S 1 1 S, 1 0, H, IJlJlEll M, P ■ In partic- 
ular, the Unruh effect [la LLZI, UM, CDl -which consists in 
the emergence of noise when an accelerated observer is 
describing Minkowski vacuum using Rindler coordinates 
-affects the possible entanglement that an accelerated ob- 
server Rob would share with an inertial observer Alice. 
The first question to be answered by theory is how much 
entanglement degrades due to Rob's acceleration a. 

Some partial answers are in the existing literature 
@, @, El • All of them share the shortcomings inherent in 
the single mode approximation (SMA) [l|, |2| , which con- 
sists in considering only one mode of sharp momentum 
in the analysis of Unruh degradation. For scalar fields, 
a Minkowskian maximally entangled state becomes sep- 
arable in the limit a — > oo, i.e. the Unruh effect com- 
pletely destroys entanglement. This is a consequence of 
the excitation of an unbounded number of modes as Rob 
accelerates. Contrary to this, finite correlations survive 
the limit a — > oo when considering fermion fields. Pauli 
exclusion principle -which bounds the maximum num- 
ber of possible excited modes- has been argued as the 
rationale for this [6(. However, that argument was only 
applied under the (somewhat unphysical) SMA. 

We will show that an unbounded number of modes 
become excited by Unruh effect even for fermion fields 
if we relax SMA, and so, the above argument ceases to 
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be plausible. Here a fundamental question arises; does 
fermionic statistics protect the entanglement? Or is this 
a mere artifact emerging from the SMA?. In this paper 
we shall show that such entanglement survival is funda- 
mentally inherent in the Fermi-Dirac statistics, and that 
it is independent of the number of modes considered, of 
the maximally entangled state we start from, and even 
of the spin of the fermion field studied. 

We will proceed step by step for the sake of clarity. 
In section [TT] we will introduce the Unruh effect and 
its impact on entanglement when one of the partners 
is accelerated. In section IIIII we will express the mul- 
timode vacuum and one particle state in the coordinates 
of an accelerated observer and for two different kinds of 
fermionic fields (A Dirac field and a "spinless" fermion 
field). After that, we will analyze entanglement degra- 
dation for two very different kinds of maximally entan- 
gled state of a Dirac field: the case of vacuum entan- 
gled with one particle state in section IIVI and the case 
of a spin and momentum maximally entangled state in 
section [Vj Both cases were considered under SMA in 
[l5| . Here, we get that, even for the radically different 
final states obtained in each case, after non-trivial com- 
putations entanglement degradation ends up being the 
same for both and, more importantly, it is independent 
of the number of modes considered. Then, in section fVTl 
we will investigate degradation for a maximally entan- 
gled state of a "spinless fermion" field (considered under 
SMA in [6|). Here, the field and occupation numbers al- 
lowed are completely different from the previous cases, 
hence, the final state is, as well, notably different from 
them. However, the dependence of entanglement on a 
turns out to be exactly the same as in the spin 1/2 cases 
analyzed before. Finally, we will trace back this behavior 
to fermionic statistics (which all the cases share). Specif- 
ically, fermionic statistics translates into a peculiar struc- 
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ture in the density matrices for Rob, which is responsible 
for those striking coincidences. 

Summarizing, although it is true that now we have 
an infinite number of excited fermionic modes, statis- 
tics counterbalances that effect and allows entanglement 
preservation even at the limit a — *■ oo, contrarily to the 
intuition we would get from the bosonic case. The mean- 
ing of this remaining entanglement and its relation with 
information content of black holes is discussed in the con- 
clusions. 



II. UNRUH EFFECT AND ENTANGLEMENT 
DECOHERENCE 

The Unruh effect appears when we try to describe fields 
in the frame of a non-inertial observer. When this hap- 
pens, an accelerated observer of the minkowskian vac- 
uum, would observe a thermal particle distribution of 
temperature Tjj = ha/2nkBC. 

To understand where this effect comes from we need to 
start from a Minkowskian frame and consider the Dirac 
field expansion in terms of the positive (particle) and the 
negative (antiparticle) energy solutions of Dirac equation 
notated tp~£ and tp^ 

s J 

Here, the subscript k notates momentum which labels the 
modes of the same energy and s = {t 7 4} is the spin label 
that indicates spin-up or spin-down along the quantiza- 
tion axis. ak, s and bk, s are respectively the annihilation 
operators for particles and antiparticles, and satisfy the 
usual anticommutation relations. 

For each mode of frequency k and spin s the positive 
and negative energy modes have the form 

^ s = _L= u ±(fc) e ±*C*-»-M) (2) 

where uf(k) is a spinor satisfying the normalization re- 
lations ±uf(k)uf,(k) — (k /m)S ss ',uf(k)uf,(k) =0. 

The modes are classified as particle or antiparticle re- 
spect to d t (Minkowski Killing vector directed to the fu- 
ture). The Minkowski vacuum state is defined by the 
tensor product of each frequency mode vacuum 

|0)=(g)|0 fe ) + |(M- (3) 

fe.fe' 

such that it is annihilated by ak, s and bk, s for all values 
of s. 

An uniformly accelerated observer viewpoint is de- 
scribed by means of the Rindler coordinates [21( . In order 
to cover the whole Minkowski space-time, two different 
set of coordinates are necessary. These sets of coordi- 
nates define two causally disconnected regions in Rindler 
space-time. If we consider that the uniform acceleration 
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Figure 1: Rindler space-time diagram: lines of constant posi- 
tion z — const, are hyperbolae and all the curves of constant 
proper time t for the accelerated observer are straight lines 
that come from the origin. An uniformly accelerated observer 
Rob travels along a hyperbola constrained to region I 

a lies on the z axis, the new Rindler coordinates (t, x, y, z) 
as a function of Minkowski coordinates (t, x , y, z) are 

at = e az sinh(ai), az = e az cosh(at), x = x, y = y (4) 

for region I, and 

at = —e az sinh(ai), az = —e az cosh(at), x = x, y = y 

for region IV. As we can see from fig. 1, although we have 
covered the whole Minkowski space-time with these sets 
of coordinates, there are two more regions labeled II and 
III. To map them we would need to switch cosh <-> sinh 
in equations {H,©. In these regions t is a spacelike co- 
ordinate and z is a timclikc coordinate. However, the so- 
lutions of Dirac equation in such regions are not required 
to discuss entanglement between Alice and an accelerated 
observer, since he would be constrained to either region 
I or IV, having no possible access to the op pos ite regions 
as they are causally disconnected [g, |g, |2l|, |22j . 

The Rindler coordinates z,t go from — oo to oo in- 
dependently in regions I and IV. It means that each 
region admits a separate quantization procedure with 
their corresponding positive and negative energy solu- 
tions 1 {^+ and {^ k v + , ^"}. 

Particles and antiparticles will be classified with re- 
spect to the future-directed timelike Killing vector in each 
region. In region I the future-directed Killing vector is 

dl=^d i+ ^d s = a(zd- t +td s ), (6) 



Throughout this work we will consider that the spin of each mode 
is in the acceleration direction and, hence, spin will not undergo 
Thomas precession due to instant Wigner rotations 0, |20| . 
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whereas in region IV the future-directed Killing vector is 

This means that solutions in region I, having time de- 
pendence ip k + ~ e ~ lk ot w ith ko > 0, represent positive 
energy solutions, whereas solutions in region IV, having 
time dependence ^ + ~ e ~ lk ot w ith ko > 0, are actually 
negative energy solutions since d\ v points to the opposite 
direction of c% @, H3| • As I and IV are causally discon- 
nected ipjX and ipuZ only have support in their own 
regions, vanishing outside them. 

Let us denote (c/,fe, s , cj fe s ) the particle annihilation 

and creation operators in region I and (dj, k , s , d\ k s ) 
the corresponding antiparticle operators. Analo- 
gously we define (civ,k,s, c\ vk s , d IV ,k : s, d\ vk s ) the par- 
ticle/ antiparticle operators in region IV. 

These operators satisfy the anticommutation relations 
{ c R,fc,s' c R',fc', s '} = $RR'hk>S ss > where the subscript R 
notates the Rindler region of the operator R = {I, IV}. 
All other anticommutators are zero. That includes the 
anticommutators between operators in different regions 
of the Rindler space-time. 

Taking this into account we can expand the Dirac field 
in Rindler coordinates analogously to (JTj) : 

^ = E [ d3k + 4,fc,X- + °iv^Z + 

s J 

+4 v ,kAy) • ( ? ) 

Equations ([T]) and (J7J) represent the expansion of the 
Dirac field in its modes in Minkowski and Rindler coordi- 
nates respectively. We can relate Minkowski and Rindler 
creation and annihilation operators by taking appropri- 
ate inner products and computing the so-called Bogoli- 
ubov coefficients @, El US [|2 : 

a fcjS = cosr cixs - e^sinr djv,-k,-« 

b Ls = cosr 4v,k,s + e~ l0 sinrc/,_ fc ,_ s (8) 

where 

tanr = e^ (9) 

and cj> is a phase factor that will turn out to be irrelevant 
for our purposes. Notice that since in Rindler regions 
I and IV the temporal Killing vectors pointing to the 
future have opposite senses, all the magnitudes that are 
not invariant under time reversal change among regions. 

It is shown in the literature 1 1 51 1 that Unruh ef- 

fect provokes decoherence of entangled states when one 
of the partners is describing the system from an accel- 
erated frame. As it can be thoroughly seen in this 
comes about because the fact of accelerating introduces 
an horizon in the space-time and, as Rob is always con- 
strained to region I or IV of the Rindler space-time, for a 
quantum description of Rob's subsystem it will be neces- 
sary to trace over the region causally disconnected from 
Rob. 



It is this partial tracing which transforms the total 
state (which would be pure for an inertial observer) into 
a mixed state whose entanglement decreases as Rob ac- 
celeration increases. This phenomenon has been called 
Unruh decoherence in the literature. We are not tackling 
here the complete problem of expressing the vacuum and 
the one particle state for the different fields beyond the 
SMA, because its calculation, although easy, may take 
a bit long and it is detailedly done in MS]. Instead we 
will start from the results obtained in [l5| for the mul- 
timode vacuum expressed as a squeezed state in Rindler 
coordinates. 

Notice that since the observer Rob is accelerated, his 
possible measurements are affected by a Doppler-like ef- 
fect. A discussion of this effect and how it would affect 
our study is given in [lj|. 



III. VACUUM AND 1-PARTICLE STATES OF 
FERMIONIC FIELDS BEYOND SMA 

In this section we shall go beyond the single mode ap- 
proximation to build the vacuum state and the 1-particle 
excited state for two very different kinds of fermionic 
fields: First a Dirac field and then a spinless fermion 
field. Both kinds of fields were analyzed under the SMA 
in previous literature (the spinless case in Q and the 
Dirac field in [l5[ whose notation we will follow in this 
paper). 

To begin with, let us consider a discrete number n of 
different modes of a Dirac field k±, . . . , k n , labeling with 
Si the spin degree of freedom of each mode, so Minkowski 
multimode vacuum should be expressed as a squeezed 
state in Rindler coordinates which is an arbitrary super- 
position of spins and momenta as it is discussed in [151 ] 

l°> = E E c ™,..,s mM ,..Mjsl;::± H z H IV 

fei,.. -,k m 

(10) 

Where, the notation is 

\i) x = l s i?^ii ■ ■ ■ > s ijki) j \ — si,— k\; . . . ;— Si,— ki) IV 

(11) 

with 

|fci,si; . . . ;fc m ,s m ) J = c\ km Sm . . . c\ MiSi \0)j (12) 

The label outside the kets notates Rindler space-time re- 
gion, and the symbol £ is if {ki, Si} = {kj, Sj} for any 
i =/= j, and it is 1 otherwise, imposing Pauli exclusion 
principle constraints on the state (quantum numbers of 
fermions cannot coincide). 

Due to the anticommutation relations of the fermionic 
operators, terms with different orderings are not inde- 
pendent. So, without loss of generality, we could choose 
not to write all the possible orderings in (TTUJ) selecting 
one of them instead. In this fashion we will write the 
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elements (jlll) with the following ordering criterion: 

h = k l+ i =>■ Si =T, =4 ■ 



(13) 



The coefficients C m are constrained because the 
Minkowski vacuum should satisfy ak ,s |0) = 0, Vfco; So- 
In [l5| we showed that imposing this constraint translates 
into 



C m = C°e im * tan" 



(14) 



Where tanr = exp (— irkoc/a). C m is independent of 
Si and fej. Therefore, we obtain the vacuum state by 
substituting (TT4")) in (TIT))) and factoring the coefficients 
out of the ki, Si summation. 



10) 



2n 

E 



C h 



ST eh 



'XV 



(15) 



fei, 



The only parameter not fixed yet is C°. We can fix it 
imposing the normalization of the Minkowski vacuum in 
Rindler coordinates (0|0) = 1, [fa, |6j, [lfa] which implies 



\C°\ =cos 2 



(16) 



Where we have taken into account our ordering choice 
explained above. This expression is explicitly derived in 
the appendix [A"l 

Eq. (fT6|) gives the value of C° except for a global phase. 
Next, the 1-particle state can be worked out translating 
the Minkowski one particle state \k,s) = a\ |0) into 
Rindler coordinates 



\k,s) 



Where 



2n-l 

E 

m=0 



A m E fcw 1^; fc ' s >/ l™>/v (i7) 



Si I 

fcl, 



= (C"» cosr + C ,m+l e -# sinr ) 



(18) 



and the notation \rh;s,k)j, consequently with pip , 
means the ordered version of |si,fci; . . . ;s n ,k n ] k, s) I . 

Another different kind of field that we are going to 
consider appears by neglecting spin keeping the fermionic 
statistics (like considering Grassman scalar fields). This 
kind of field is used under the SMA in the literature Q. 
Here we will relax such approximation and will analyze 
Unruh dccohcrence when we allow n different momenta 
ki. Barring spin, the Minkowski multimode vacuum state 
would be expressed as 



|0) 



E ^ *. 

m=0 k\,...,k m 



I IV 



(19) 



where, in this occasion |m) 7 \rh) 



XV 



h 

i "'mil 



Using the same procedures as for the 



spin 1/2 case (|T0|) we can prove that all the coefficients 
are independent of fc, and can be related to C° as in (| 14[) , 
C m = C e im( t> tan m r. We can now fix C° imposing the 
normalization relation (0 |0) = 1 giving 



C° 



,m=0 



-1/2 



where 



X, 



= E ^i,...,k m = 



fcl , • --jfem 

Such that ([20]) can be simplified to 



(20) 



(21) 



C° 



E 

.m=0 



tan 2 " 1 r 



-1/2 



cos" r (22) 



Finally, the one particle state a\ |0) is 



\k) = ]T i m E &i.-.*».k l™> k )i \™)iv ( 23 ) 

m=0 k\,...,k m 

where A m has the expression IT8|) but substituting C m 
by C"™ 



IV. ENTANGLEMENT DEGRADATION FOR A 
|00> + 11) ENTANGLED STATE OF A DIRAC 
FIELD BEYOND SMA 

In the following we will analyze Unruh entanglement 
degradation in various settings corresponding to different 
maximally entangled states of fermion fields. First we 
consider the following state in Minkowskian coordinates 



1 



V2 



(\0) A \0) R + \k A ,s A ) A \k R ,s R ) R ) (24) 



The density matrix for the accelerated observer Rob is 
obtained after expressing Rob's state in Rindler coordi- 
nates -which means using (fT0|) and (fl7|) in Rob's part of 
(IM)) - and afterwards, tracing over Rindler's region IV 
since Rob is causally disconnected from it and he is not 
to extract any information from beyond the horizon. Fol- 
lowing this procedure we obtain the density matrix p = 



MEK E &±|0>J™>z<0U<" 

m=0 S!,...,s m 
fcl ,. . . ,k m 

2n-l 

+ E {p? E fctiri^Ni^^ 

m=0 si s m 

k± . . . .,fc m 

2n-l 

x (m; k R , s r \^+J2(d? ]T I**, s A ) A 

m=0 si,...,s m 
fcl ,...,km 

x |m; fefl, s R ) / (k A , s A \ A (m; s R |j J + (H.c.) non _(25) 

' - diag. 
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Where (H.c.) non -diag. means Hermitian conjugate of only 
the non-diagonal terms and 



, tan 



(26) 



with i = 0, 1, 2. The derivation of (|25[) can be found in 
the appendix [Bj 

Notice that as Rob accelerates, the state becomes 
mixed with all the available modes (hi, . . . , k n ) excited. 
This contrasts with the Minkowskian state where 
only one Rob mode is excited (k R , sr). Notice also that 
the Hilbert space dimension has changed from 2 x 2 to 
2n x 2n. 

We will compute the negativity as a function of a as a 
measure of the state entanglement, the negativity is the 
sum of all the negative eigenvalues of the partial trans- 
pose of p. 

The partial transpose of (f2"5| has a 2 x 2 and lxl 
blocks structure. Each eigenvalue in the lxl blocks 
is non-negative (since Df 1 > 0), so we are interested in 
the 2x2 which are the ones that may have negative 
eigenvalues. These 2x2 blocks expressed in the basis 

r 2n-l 

\ |0) A |m; k R , s R )j , \s A , k A ) A \fh)j \ (27) 

1^ ) m— 



are of the form 



m+l 








(28) 



There is no element because it goes with 

\k A ,s A ) A \m;k Rl s R ) I (k Al s A \ A (m;k R ,s R \ I which can- 
not have any element within our block as Pauli exclusion 
principle imposes k R , s R £ {k % , s % } i=1 ^ m . 

Each 2x2 block of ([28| appears a number of times 
B m given by 



B m = 



2n- 1 



(29) 



The derivation of this formula can be found in the ap- 
pendix El 

Using ([26]) . the negative eigenvalue of each block can 
be expressed 



IA-I = l\C°\ 2 tan 2 



(30) 



where Cq is given by (|16p . Therefore, the negativity is 
expressed as the summation of the negative eigenvalue 
of each block |A~| multiplied by the number of times 
B m that that block appears in the partially transposed 
density matrix. The summation of the series is 



M = ]T B m \\; 



m=0 



m=0 \ 



(31) 



but this result can be easily simplified to 



TV 



1 



cos 2 r 



(32) 



which is independent of the number of modes n that we 
have considered. This surprising result shows that, even 
if we consider more than one mode in our field analysis, 
the entanglement degradation due to Unruh effect is the 
same as considering only one mode as it is done in [l5j |. 
In other words, despite the fact that all the available 
modes are excited when Rob accelerates (|25|) . the quan- 
tum correlations behave as if we were considering only 
one possible mode for the field. This is a consequence of 
the peculiar structure of the density matrix for Rob, be- 
ing the fermionic nature of the field the final responsible 
of this structure (f2"T)h 



V. ENTANGLEMENT DEGRADATION FOR A 
SPIN AND MOMENTUM ENTANGLED STATE 
OF A DIRAC FIELD BEYOND SMA 

If instead of (|24p we start from a Bell momentum-spin 
state in Minkowski coordinates 



I*) 



k A , S A ) A \k R , s R ) R + \ k A ,s\) A \ k%, s R ) R ) 

(33) 

As it can be seen in the appendix [B] the density matrix 
for Rob takes the form p = 



2n-l 



71=0 Sl,...,S m 

k\ ,...,k m 

I ~ ; 1 II s-ki,...,k rn ,k'^ t It 2 2\ I ~ I 2 2\ 

x \ m j s Ii\i + ? Sll ...,a m ,g2 s A/ A | to ; fc fi> S R ) 



Ma 



(kl,s 2 



A i S A\ A 



(m; k%, s 2 



x |)7i; fc^j, Sfl) j (k A , s 



R>°R\l ■ * 8l 
A; a A\ A (" l > ^R' a R\l 



[Tfl\ kr>, Si 



(H.c.) non _ (34) 

diag. 



Analogously to (|23)) . the partial transpose of ([3^)1 has a 
2x2 and lxl blocks structure. Again, we are interested 
in the 2x2 blocks -the ones that may have negative 
eigenvalues.- These blocks expressed in the basis 



I ~ 7 2 2 

m k 



are of the form 



2 2 ill 2n ~ 2 

• " °r)t ' \ S A' k A ) A \rh, k R , s R ) > 

J m— 

(35) 



±Dl 



±D r 9 







(36) 



Notice that there is no diagonal elements in the block 
because the terms that would go in the diagonal are for- 
bidden by Pauli exclusion principle, which imposes that 

This time, each 2x2 



°R> ^_R' S R & {fci, Si} i=1 
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block of the form (|36|) appears a number of times B' m 
given by 



B' = 



In- 2 



(37) 



(See appendix |Aj. The negative eigenvalue of each block 
is 



A™ — 



D™ cos 4 ™- 2 r 



tan 2m r (38) 



would have the form 



1 / & n+1 ±D\ n 

2 v ±bf 



(44) 



The main difference with (|28[) is that Co a different value 
(I2"2"j) instead of C° given by (H]). Here, I)™ does not 
appear because Pauli exclusion principle imposes that 
k R ^ {ki} i= i m - Now, each 2x2 block multiplicity is 



where Co has been substituted by (fT6|) . Therefore, the 
negativity results 



M=J2 B 'ra\K 



cos 4 "" 2 r 2 ^ 2 (In - 2 



E 



2 ^ \ m 

m=0 



tan 2m r 



(39) 



(40) 



This can be readily simplified to 

Af = - cos 2 r 
2 

Strikingly we run into the same simple result as above 
(|32|) . Even starting from a spin Bell state, the entangle- 
ment is degraded by Unruh effect in the same way as in 
the previous case. 



VI. ENTANGLEMENT DEGRADATION FOR 
|00) + |11) ENTANGLED STATE OF A SPINLESS 
FERMION FIELD BEYOND SMA 

Now we can go one step further neglecting spin and 
consider a spinless field on which we have imposed the 
fcrmionic statistics. The maximally entangled state of 
the vacuum and one particle in this setting 



^) = ^(\0) A \0) R + \kA) A \k R ) 1 



(41) 



As it is discussed in appendix [Bl this leads to the fol- 
lowing density matrix for the accelerated observer Rob 
after using expressions (fTi?)) and (f2"3"|) and after tracing 
over Rindler's region IV p = 

n n— 1 

m—0 ki,...,k m m—0 

Zk 1 ,...,k m ,k R \0} A \m} I (k A \ A (rh;k R \ I + D™ £ 



£,k u ...,k m ,k R \k A ) A \m; k R )j (k A \ A (m; k R \j 



ki ,. . ■ : k 7n 

(H.c.) non _ 

diag. 

(42) 



where D\ is given by the expression (f2T>)) but substituting 
C° by C . 

Analogously to (|23|) and (13~4"|) , The partial transpose of 
(g2) has a 2 x 2 and lxl blocks structure. The 2x2 
blocks expressed in the basis 

\\Q) A \ih ] k R ) I ^k A ) A \m) I \ (43) 



B'' 



n- 1 



(45) 



(See appendix). The negative eigenvalue of each block is 
given by the same expression (I30[) but Co is now given 
by (|2"2"|) , which is to say 

| A- | = 1 1 C° | 2 tan 2 " r = i cos 2m tan 2m r (46) 



We can compute the negativity 



n — 1 n — 1 

=]>>;;jA m | = - cos 2 " r]T 



m=0 



n — 1 



tan 2m r 



(47) 

At this point, the reader might not be surprised by the 
resulting negativity after straightforward simplification 



N = — cos 2 r 
2 



(48) 



which is the same result as in the cases ([M]) and (|33p . 
Again, entanglement degradation due to Unruh effect is 
the same as considering one mode of a Dirac field fl5T ]. 



VII. CONCLUSIONS AND COMMENTS 

Let us summarize our results so far. We have stud- 
ied entanglement degradation by Unruh effect as Rob 
accelerates beyond the single mode approximation and 
three different Minkowskian maximally entangled states: 
1) Vacuum- vacuum plus one-particle-one-particle maxi- 
mally entangled state of a Dirac field, 2) spin-momentum 
Bell state for a Dirac field, 3) Vacuum-vacuum plus one- 
particle-one-particle maximally entangled state of a spin- 
less fermion field. In spite of the essential differences 
among these states, the negativity degrades in exactly 
the same way for any acceleration. This result may look 
surprising considering that this is the same degradation 
obtained under the single mode approximation 0, [l5j 
but as it is discussed in this paper, this is an outcome of 
fermionic statistics. 

In the bosonic case acceleration excites an infinite num- 
ber of modes, and this completely degrades the entangle- 
ment in the limit a — * oo. Although one could expect 
the same behavior here -as an infinite number of modes 
is also excited when we let n — > oo- our results show 
that some entanglement is preserved for a — > oo. It is 
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remarkable that the entanglement degradation coincides 
for all the different cases considered, with independence 
of the number of modes n. 

This striking result can be traced back to the fanciful 
block structure of Rob density matrix, which produces 
the same negativity even when the characteristics of the 
entangled states (and even the field) change. The culprit 
of this structure is fermionic statistics, (as we have dis- 
cussed after ([251) . ([55]) . ijS])) which is responsible for the 
identical, and somewhat unforeseen, negativity behavior. 
This is a global feature of maximally entangled states of 
fermion fields and not a consequence of the specific cases 
chosen and the number of modes considered. 

So, Af — > 1/4 when a — > oo, and this happens indepen- 
dently of the number of modes of the field that we are 
considering, of the starting maximally entangled state, 
and even of the spin of the field which we study. What 
all the cases have in common is the fermionic statistics 
itself, so, widening the margin for Unruh degradation for 
fermionic fields will not affect entanglement degradation. 

Notice that a very different scenario would come from 
a setting in which we erase partial information for the 
state as Rob accelerates (e.g. angular momentum). In 
that case, it was shown that entanglement degradation 
is greater than in the cases where all the information is 
taken into account [l5| , but this has more to do with this 
erasure of information than which the fermionic nature 
of the states. 

One question immediately arises from these results; 
Are the remaining correlations purely statistical? In 
other words, does any useful information really survive 
the limit a — ► oo?. As all the states undergo the same 
degradation, everything points that statistics is the only 
information which survives this limit. 

Furthermore, the limit a — * oo can be understood as 
considering an observer moving in a trajectory arbitrarily 
close to the event horizon of a Schwarzschild black hole 
Q. So, even if Alice is free falling into a black hole and 
Rob stands at the event horizon, a fixed degree of entan- 
glement survives to Unruh decoherence. Apart from the 
interest of describing the entanglement between accel- 
erated observers, the regularity and universality of our 
result (N — (1/2) cos 2 r) could be a useful hint in the 
solution of the information paradox in black holes and 
deserves further investigation in future works. 



Appendix A: DERIVATION OF Co AND THE 
COMBINATORY FORMULAE 



To derive Co except for a global phase, we impose the 
normalization of the vacuum state in Rindler coordinates 
(0 |0) = 1, from (|15p. we see that this means that 



where 



T m tan 2m r 



J2 T 2n _ m tan 2m r 



.m=0 



m— n+1 



J-TO — ? si,...,s m 
si ,. . . : s m 
ki , . . . : k m 

Now, we are going to show that (|A2j> , has the form 



-1/2 

(Al) 
(A2) 



p=0 



n — p 



m - 



2p. 



ym. — 2p 



(A3) 



To see how this expression comes from Pauli exclusion 
principle, we have to read p as an index that repre- 
sents the number of possible spin pairs (ki = fcj-j-i, s,: =| 
, Si+i =1) which can be formed, and goes from to the 
integer part of to/2, and then 

• The combinatory number (j^Zf ) represents the 
possible combinations of modes that can be formed 
taking into account that p different momenta ki are 
not available since they are already occupied by the 
p pairs. Hence, it is given by the combinations of 
the n — p available momenta taken m — 2p at time, 
since m—2p is the number of free momentum 'slots' 
(the total number of different momenta m minus 
the number of positions taken by pairs 2p). 

• The combinatory factor (™J represents the different 
possible combinations for the configuration of the p 
pairs, which have n possible different momenta to 
be combined among them without repetition and 
in a particular order. 

• The factor 2 m_2p represents the possible combina- 
tion for the spin degree of freedom of each mode. As 
a spin pair only admits one spin configuration, only 
the unpaired modes will give different spin contri- 
butions, so the factor is (2S + l) m ~ 2 P giving the 
formula (jA"2"j) 

After some lengthy but elementary algebra we can see 
that 
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T — 



(A4) 



and using the property ( " .) = (?) , we can express (|Al 
as 



E ( 2n ) ' 

,m=0 ^ 



-1/2 



cos 2 ™ r 



(A5) 
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which is equation ([16]) . 

Now we will do the same for the equation (|29[) . This 
formula takes into account the number of two by two 
blocks of the form ([25]) . Taking a look at the basis in 
which those blocks are expressed (|27|) . we can see that 
the expression for B m is given by two terms: 

• The number of possible combinations of to modes 
with n possible different momenta ki and two pos- 
sible spins Si according to Pauli exclusion principle 
as in (EH). 



A negative contribution which comes from exclud- 
ing those combinations in which {kR, sr} coincides 
with any {ki, Si}, which means excluding the num- 
ber of combinations in (|A3|) which have one of their 
values fixed to {ki, Si} = {kR, sr}. This number is 
given by the combinatory number ( n ~i) provided 
that m > and it is zero if m = 0. 



To see where this negative contribution comes from let 
us assume that it is {ki,st} the mode which coincides 
with {kft, sr} we will have 2n — 1 possible choices for 
each {kj^i, Sj^i} (2 values for s and n for k excepting 
ki , Si due to Pauli exclusion principle) . This happens for 
all the combinations of all the possible values {kj,Sj} 
with j ^ i. Hence, as there are m modes and one of 
them is fixed ki = kn,Si — sr, we have to consider the 
combinations of 2n — 1 elements taken m — 1 at time. 

If to > n the situation is equivalent to having m! — 
2n — to. Since having more modes to than possible values 
of ki we are forced to have n—m pairs and we lose freedom 
to combine the available modes. 

Now if we compute 



B rl 



T 



2n 

TO - 



2n- 

TO - 



After some basic algebra we obtain 

"In - 1 

TO 



B,, 



(A6) 



(A7) 



which is expression 

The derivation of expression ([57)) is quite straightfor- 
ward considering the one above. Looking at the basis 
of the 2x2 blocks (l35j) we can see that this case would 
be exactly the same as the previous one but now {ki, Si} 
cannot coincide neither with {k R ,s R } nor {k R , s R }. Re- 
peating the same reasoning as before we have to do three 
operations as follows 

• Discounting the combinations which have a coin- 
cidence {ki,Si} — {k R , s R } from the total number 
(IA3|) and obtain the expression (|A7|) 

• Subtracting the combinations with coincidences 
{kj,Sj} — {k R , s R } 

• Taking into account that we have subtracted twice 
the cases in which we have double coincidences, we 
need to add the number of double coincidences once 
to compensate it. 



The number of cases with double coincidences (which 
require to > 1) is the combinatory number ( ™~ 2 ), as 
we have 2n possible spins and momenta minus the two 
fixed possibilities ({ki,Si} — {k R ,s R } and {kj,Sj} — 
{k R , s R }) and to modes being 2 of them fixed. Taking 
this into account 



(2n- 


; 




/2n- 


0- 


,2n- 


I) 


V TO — 






V TO — 




V TO — 





This expression can be simplified to 



(2n- 







(2n - 2\ 


V TO — 






\ m J 



(A8) 



(A9) 



which is (|57j) 

For the spinless fermion field, equation (|21|) . which 
have the form 



E 

k± , . ...kr, 



(A10) 



corresponding to the possible combinations of m values of 
ki imposing that ki ^ kj if j =/= i (which is the translation 
of Pauli exclusion principle to spinless modes). This can 
be readily obtained taking into account that the n pos- 
sible values of ki should be combined without repetition 
in a particular ordering of the to modes, so the possible 
combinations are simply the combinatory number (^) 

Equation (|45|) can be easily obtained taking into ac- 
count that the number of 2 x 2 blocks (|43[) is given by 
the number of mode combinations allowed by Pauli prin- 
ciple (|A10[) . subtracting the terms having ki = kR. The 
number of possible kj values allowed for the rest to — 1 
modes having fixed ki = kR is n — 1, so the number of 
combinations we must subtract is the combinatory num- 



ber 



obtaining 



B^ = 



which is (US' 



n- 1 

TO — 1 



n - 1 

TO 



(All) 



Appendix B: DENSITY MATRIX 
CONSTRUCTION 



In this appendix we will derive expressions ([25]) . (|34l) . 

for the density matrix of the system Alice-Rob. 
Using expression (fT5|) we see that the Alice-Rob 
Minkowskian operator Poo = |0;0)(0;0| when Rob is ac- 
celerating translates into 

2n 2n 

poo - EE r ( c 'r E &± E4'.'.:? 

ki : ....k l 

)a\™)i\™) iv (1\iv(1\i(0\ a (Bl) 



m=0 1=0 
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where 

= (k' 1 ,s[;...]k , l ,s'i\ I = {Q\ I c Iyk ' i y i ...c Iikli 
{—k ll —Si, . . . ; — ki, —Si\j 



IV 



\IV c /V,-fei,-si ■ ■ • ClV.-k'^ 



(B2) 



In expression (jBll) , and below in (|B3|) . bras and kets 
refer to Alice's mode in Minkowski coordinates and Rob's 
mode in Rindlcr coordinates. 

Now, using expression (JTTJ) we can write the operator 

P %] — \h l a' • lc* a' \lh3 a- 7 ■ P <J I 
r ll — \ A ' A ' H> i?7 \ A ' A ' ' R\ 



p'J 



2n-l 2n-l 

E E 

m=0 Z=0 



Sl , . ..,Sr, 

ki , . . .,/c„ 



7V 



Si,-,*, 

/ci ,...,k'. 



(B3) 



where A™ is given by (jT8j) . 

Notice that the objects m; k R , s^) 7 represent the ap- 
propriate ordering of the elements inside with its sign, 
taking the criterion (|13[) into account. 

Now we can use expressions (fl~5)) and (IT7)) to obtain 
the operator Poi = |00)(fc J 4, sa', k R , s R \ as it is expressed 
when Rob is describing the world using Rindler coordi- 
nates. 



2n 2n-l 



m=0 Z=0 



Si,...,S 7 



s 1 ,...,s l 
ki ,...,k[ 



x I™)/ ( l Uv(l; kn, S R \! (0\ a 



(B4) 



After obtaining the expressions for the operators 
-Poo, Pu, Pqi we can write the density matrix associated 
with the state (1241) in Rindler coordinates for Rob, 



(B5) 



P=\ (Poo + Poi + P i + Pn 



Where for Pl\ we are considering {k R , s R } — {k R , s R } = 
{k R , s R } and {k\, s\} = {k 3 A , s 3 A } = {k A , s A }. 

We can do the same to obtain the density matrix as- 
sociated with (|3"31 in Rindler coordinates for Rob 



P=^(PX+P%+Pl? + (P^) (B6) 

Now, we must consider that, as Rob is causally discon- 
nected from Ridler's region IV, we should trace over that 
region to obtain Rob's density matrix. Hence, we need 
to compute the trace over IV for each of the previous 
operators (jBTj) . (E3J) . |B4|) . 

Taking this trace is actually quite straightforward tak- 
ing into account the orthonormality of our basis once we 
have chosen one particular ordering criterion (|13|) . 

(fh\rh') IV = 5 mm > (Ss^Sk^ . . ■S Sm , s ' m Sk m ,k' m ) (B7) 



Hence, 



Tr/yP o= E (™'\iv P oo \m') 



iv ■ 



(B8) 



m'—0 



Using (|B7[) only the diagonal elements in region TV sur- 
vive and (IB8I) turns out to be 



Tr IV Poo = £ \C m \ 2 ]T ^:;;± \0) A \m) I {m\i 



m=0 Sx,.,.,s n 
k\ ,. . 



(B9) 

which, substituting C m as a function of C° using |T 
and then (f2"6| , is expressed as 

2n 

Tr IV Poo = £ D™ &± I°)a (m|, 



m=0 si,...,s 7T 
i , . . . , k „ 



(BIO) 



Now we will compute the trace 

2n 

Tr/y PH = £ (m'U P H \™')iv ( B11 ) 



m'=0 



2n-l 



Si ,. . . ,5^ 



x \k A ,s A ) A \rh;k R ,s R )j (rh;k R ,s R \i (k A ,s A \ A 

(B12) 

substituting C m as a function of C° (combining (|18p and 
([H])) we can express |A m | 2 as 



o, / sin" r \ ,„ l0 tan 2m r 
|C | 2 tan 2m r cosr+— — ) = |C | 2 — = D 



COS T 



Such that 

TrivPii = E D 2 Y,tsl;::±;s«\k A ,SA) J 



2 

(B13) 



2n-l 



=0 s 1 ,...,s„ 
fci 



x |m; 4, (m; 4, s^|j (k A , s\ jBU) 

When {fe^^jj} = {k R ,s R } = {k R ,s R }, {k\,s A } = 
{k A ,s 3 A } = {k A ,s A }. 

However, in the general case i ^ j it would be 



Tr 



pii = E^ n E C 



m^O si,.. M s T 
fci ,. . . ,k r 



'■ ><-'**•»"*>)* \k A ,a A ) 



A 



x \rh; k R , s R )j (m; /c fl , s_r|/ (fc A , s^U (B15) 
Now, let us compute the trace 



2» 



Tr/vPti = ^2 (™'\iv p ai I" 1 ') 



i\- 



(B16) 
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2n-l 

Tr /y Poi = E E l™>/ \™)iv 

fei , . ..,k m 

x(l\MhkR,8 R \i(0\ A (B17) 

from HU) and (O we see that the product C m (A m )* is 
real and has the expression C m (A m )* = 

|C | 2 tan 2m r fcosr+— ^1 = |C | 2 ^^ = D? 
\ cos r / cos r 

(B18) 

so that 

2n-l 

Tr /y Poi = E D ? E |0>a l™>/ \™)iv 

m=0 si,...,s m 
fcl ,...,fe m 

X^l/y^fcfl, fl fl |/<Q| A (B19) 

Now we can compute Rob's density matrices for each 
case tracing over IV in expressions (|B5[) and (|B6p . First 
the matrix (|B5[) is, after tracing over IV, 

Tr/y p = - (Tt iv Poo + Tr IV P 01 + Tr IV P^ + Tr IV P») 

(B20) 



substituting expressions (|B10|) . (|B15|) . (|B19|) we get ex- 
pression (J2U). 

Now, concerning (|B6j) 



Tr IV p=- Trjy (P lT + P 22 + P\l + (P^^t) (B21) 



Substituting expressions (|B14|1 and (|B15[) we obtain ex- 
pression (|34|) . 

The derivation of (|42| is completely analogous to (I25[) , 
taking now into account that we have C m and D m in- 
stead of C m and D m and that we have no spin degree of 
freedom. Notice that, even though the structure of P2"| is 
completely analogous to the structure of (|25[) . and there- 
fore, repeating the derivation will add nothing to this 
appendix, these density matrices are completely differ- 
ent due to the different dimensions, the different values 
of C° and C° and the number of 2 x 2 blocks which give 
negative eigenvalues. 
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